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We calculate the chameleon field potential for ultracold neutrons, bouncing on top of one or 
between two neutron mirrors in the gravitational field of the Earth. For the resulting non-linear 
equations of motion we give approximate analytical solutions and compare them with exact nu- 
merical ones for which we propose the analytical fit. The obtained solutions may be used for the 
quantitative analysis of contributions of a chameleon field to the transition frequencies of quantum 
states of ultra-cold neutrons bound in the gravitational field of the Earth. 
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I. INTRODUCTION 

A chameleon field has been suggested to drive the cur- 
rent phase of cosmic acceleration for a large class of scalar 
potentials. The properties depend on the density of a 
matter to which it is immersed. Because of this sensitiv- 
ity on the environment, the field was called chameleon, 
and it can couple directly to baryons with gravitational 
strength on Earth but it would be essentially massless on 
solar system scales [U |2] • 

An interaction of a chameleon field to an environment 
with a mass density p and its self-interaction are de- 
scribed by the effective potential V e e(4>) [Sj 



V eS (<j>) = V^)+pe^ M , 



(1) 



where f3 is a coupling constant and M = 1/ \/8ttGn = 
2.435 x 10 18 GeV is the Planck mass. 

As has been pointed out by Ref. [3J, ultra-cold neu- 
trons, bouncing in the gravitational field of the Earth 
above a mirror can be a good laboratory for testing the 
existence of a chameleon field. There has been found 
a solution of the equations of motion and limits for the 
coupling constant /3 can estimated from the contribu- 
tion of a chameleon field to the transition frequencies of 
the quantum gravitational states of ultra-cold neutrons 
bouncing in the gravitational field of the Earth. Such 
resonant transitions between quantum states of a neu- 
tron in the gravity potential have been measured by the 
qBounce Collaboration [5]. In this experiment, ultra-cold 
neutrons form gravitationally bound quantum states be- 
tween a neutron mirror on bottom and another neutron 
mirror on top with a relative distance d. The upper neu- 
tron mirror has a rough surface and thereby filters out 
higher, unwanted states. The whole system may be vi- 
brated in vertical direction with frequency lo and ampli- 
tude a. Behind this system, the neutron transmission 



is measured. It shows a significant reduction in count 
rate, if the vibration frequency corresponds to the energy 
difference between two eigenstates. With this method, 
the transition frequencies have been determined, which 
tests Newton's Law at short distances and is therefore 
sensitive to any deviation like chameleon fields. So far, 
such chameleon fields have only been calculated in the 
infinitely large spatial region above one mirror [3j. For 
the experiment presented in [5] , it is necessary to use so- 
lutions of a chameleon field localized between two mirrors 
for the estimate of the coupling constant /3. 

As we show below the appearance of an addition mirror 
complicates the problem of the calculation of a chameleon 
field substantially. As a result we propose the analysis of 
all possible solutions of the problem under consideration. 

The paper is organized as follows: In section [TTJ the 
chameleon field for ultra-cold neutrons, bouncing in vac- 
uum in the gravitational field of the Earth above a hor- 
izontal mirror, is calculated. We show that the solution 



obtained in [3] is exact for this problem. In sections III 



the calculation of the chameleon field is adapted to the 
experimental setup used in [5]. For this purpose, we con- 
sider a symmetric geometry, where the horizontal mirrors 
occupy the spatial regions z 2 > 4-, whereas ultra-cold 
neutrons bounce in the spatial region z 2 < \. We ob- 
tain the solutions for a ch ame leon field in the spatial 
d 2 



region z < 



In section 



IV 



we solve the problem of 



a chameleon field between two mirrors numerically with- 
out any approximation and give an analytical fit of these 
solutions for arbitrary power n. 



II. THE CHAMELEON FIELD COUPLED TO 
ULTRA-COLD NEUTRONS BOUNCING ABOVE 
A MIRROR 

According to Ref.[3j[4], the potential V{<f>) takes the 
form: 
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V{<t>) - A 4 + 



A 4+n 



(2) 



2 



The scale A is chosen to be equal to A = 2.4(1) x 
10~ 12 GeV [3J. Following [3] and keeping only the lin- 
ear term in the expansion of the exponential e^^l M in 
powers of </> one can show that the effective potential has 
a minimum 



Vl s {4> wiu ) = V'{(t> w in) + P~ = Q, 

where (<f>) is a derivative with respect to (j>, at 



= A 



nA/A 3 \V(»+i) 



/"nMA"\ 



(3) 



(4) 



Quanta of a chameleon field are massive with mass de- 
fined by the second derivative of the potential V(<fi) [3J. 
A chameleon field satisfies an equation of motion [3J: 



□ = -Kff (0) = -V'(4>) + V'tenin). 



(5) 



As analysed by Brax et al. 3], a neutron mirror with 
mass density p m being parallel to the (x, y)-plane at z = 
0, occupies a spatial region z < 0, whereas ultra-cold 
neutrons bounce above in the gravity field in the spatial 
region z > with vacuum density p v ~ 0. 
For such a geometry a chameleon field depends on the 
spatial variable z only and the equation of motion takes 
the form 



d 2 ^ 
dz 2 



= V'{4>) - V'fomin), 



(6) 



where m ; n = cj> m and m ; n = <ftv for spatial regions z < 
and z > 0, respectively. In our analysis of a chameleon 
field for a system mirror-vacuum we follow [3J, but - in 
view of the following sections - locate the mirror at z = 
— ~. In this case the mirror occupies the spatial region 
z < — | > whereas vacuum fills up the region z > — | . 

After the first integration the equation of motion 
Eq. (JsJ) reduces to the form 



i /#y 

2\dz) 



C + V(4>)-<I>V' (trata), 



(7) 



where C is a constant of integration. This constant we 
obtain from the asymptotic conditions. This constant we 
obtain from the asymptotic conditions assuming that a 
chameleon field tends to its minimum and a vanishing 
derivative at z — ¥ ±oo |3J: 



lim i 

;— >-±oo 



These conditions result in 



lim ^1=0. (8) 

^±oo dz 



c = -v(4> min ) + ^ mln y'(^ min ). (9) 

Substituting Eq.Q into Eq.Q we obtain 
1 /dd\ 2 

2\£) = ^)- y (^)-^-^) y, (^)- ( 10 ) 



Thus, in spatial regions z < — | and z > 
field should obey the equations 



a chameleon 



1 /d(j)\ 2 
2\dz 



V(<t>) - v(<t>, 
A 4 



(4> - 4> m )v'{(t> m ) 

(f) n ■ 



and 



1 /d(j> 
2\dz 



+ (11) 



A 4+n ( , ,6 n 6 n+1 \ , . 



respectively, where we have used the exact shape of the 
potential V(4>), given by Eq.Q, and 4> m and 4> v are given 
by Eq. Q for p = p m and p — p v , respectively. 

The solutions of Eq.pT]) and Eq.(12| should obey the 
boundary conditions 



4>(z) 

d<j>{z) 
dz 



2-> — « — 



0(z) 

d<j)(z) 
dz 



>-%+ 



(13) 



where 0d = <fi(—d/2). Since ultra-cold neutrons are in 
the spatial region z > — |, we have to search a solution 



of Eq.( 12 ) only. As regards Eq.( 11 1, it may be used only 



for the boundary conditions. 

As a vacuum density p v ~ leads to 4> v — > oo, we take 
the r.h.s. of Eq.(12| in the limit <j) v — > oo. This reduces 

2 A 4+ 



Eq.(12) to the form 



dz 



(14) 



Since a chameleon field grows with z — > oo, in the spatial 
region z > — ~ a derivative of a chameleon field is always 
positive. This gives a differential equation 

, n , , 

72 Adz = (|) 



A 



(15) 



The solution of this equation is given by the integral 
' d \ f f <t>\2 d(f> 



V2A 

Integrating over <; 
with respect to 



AJ A 



(16) 



<f> and solving the obtained expression 
(z) we obtain 



where 1/A d = Ay/n(n + l)(A/ (j) d ) n + 2 . At z = -| from 
a boundary condition on the derivatives of a chameleon 
field 



A 4+r 



(l-(n + l) 



A 4+, 



(18) 



3 



we define <pd in terms of < 



Od = 



(19) 



III. THE CHAMELEON FIELD COUPLED TO 
ULTRA-COLD NEUTRONS BOUNCING 
BETWEEN TWO MIRRORS 



This relation is valid for arbitrary value of the coupling 
constant (3. Since the solution Eq.(17) obeys the asymp- 



totic conditions Eq. (JSj) at z — > oo and the boundary con- 
ditions Eq.( 13 1, it defines a chameleon field in the spatial 
region z > — |. 

It is obvious that if a mirror with a density p rn and 
vacuum with a density p v ~ occupy the spatial regions 
z > +| and z < +|, respectively, the solution of a 
chameleon field in the spatial region z < +| is 



Following the geometry used in [4] for the calculation 
of a contribution of a chameleon field to the Casimir force 
induced between two parallel plates perpendicular to the 
z-axis, we let ultra-cold neutrons bounce in the region 
z 2 < ^ with a density p v . The regions z 2 > ^ are 
occupied by neutron mirrors with a density p m such as 

Pm > Pv 

For such a geometry a chameleon field depends on the 
variable z only. It is described by the equation of motion 



0(«) = 0d(] 



n + 2 



1 



y/2n(n + 1) A d 



(20) 



n + 1 



where ^ = </>(+d/2) 

The solutions Eq.(16) and Eq.(20) agree well with the 
solution, obtained in |3J. Of course, the solution for a 
chameleon field in the spatial region above a mirror, i.e. 
2 > — |, is a meaningful from the point of view of the 
contribution of a chameleon field to the transition fre- 
quencies of the quantum gravitational states of ultra-cold 
neutrons [3J. Following [3] we define a contribution of a 
chameleon field to the gravitational potential of the Earth 
coupled to ultra-cold neutrons above a mirror as 



$(z) = rag \ - + ; 



M 



(21) 



where m and <? are the neutron mass and the gravita- 
tional acceleration, respectively. In comparison with [3] 
the mirror is shifts down to z = — | . 

The calculation of the contribution of a chameleon 
field, for example, to the transition frequencies of the 
quantum gravitational states of ultra-cold neutrons is re- 
lated to the calculation of the matrix elements (p\<f>(z)), 
where \p) is a low-lying quantum state of ultra-cold neu- 
trons in the gravitational field of the Earth [5J [7[ . Due 
the Airy functions, describing the gravitational states 
of ultra-cold neutrons, the main contribution to ma- 
trix elements (p\4>(z)\p) comes from the region around 
(| + z) ~ £ = (2ra 2 g)~ 1 / 3 = 5.9 pm, which is a natu- 
ral scale for neutron quantum states in the gravitational 
field of the Earth [6j [7] • In this case in the strong cou- 
pling limit /3 > 10 5 due to a relation Eq.(19) and the 
experimental density of a mirror p m = 2.51 g/cm 3 [5] 
one may find that the contribution of the first term in 



the parentheses of Eqs.(17) and (20) can be neglected 



and the solutions for a chameleon field take the form 



(22) 



independent of the coupling constant /3. This agrees well 
with the result, obtained in [3J. 



dz'< 



= V'{</>) - V (</» min ). 



(23) 



This equation is valid for a chameleon field in three spa- 



tial regions z > 



4 and z 2 < ^ 



At z = ± % a chameleon 



field has to satisfy the standard boundary conditions 



4>{z) 

d<j>(z) 
dz 



7 _1_ „ 



d<j)(z) 
dz 



(24) 



In the spatial regions z 2 > ^- after the first integration 
Eq.(23) reduces to the first order differential equation 



(see Eq.(lll) 



1 /#y 
2\dz) 



V(4>) - V{<l>„ 



(</> - 4>m) v'{4> n 

Ln+1 , 



^r(}- {n+i) k +n yy (25) 



where the integration constants are defined from the 
asymptotic conditions at z — > ±oo. Below we use Eq.(25) 



only for the boundary conditions as we have done in the 
case of one mirror in section ITU 

Before the integration of the equation of motion for a 



chameleon field in the spatial region 



< 



we have 



to accept that a chameleon field has a minimum at <p v , 
which is caused by the properties of the effective poten- 
tial V e ff((f>) but not a spatial region of a location of a 
chameleon field. Since the region of a localisation of a 
chameleon field z 2 < 4- is finite, such a minimum can be 
never reached inside. 

After the first integration the equation of motion for 
a chameleon field in z 2 < j reduces to the first order 
differential equation 



dtp 
dz 



C + V(<j>)-4>V'(<l> v ), 



(26) 



where C is a constant of integration. Since the asymp- 
totic regions z — ¥ ±oo are not reachable for a chameleon 
field in the spatial region z 2 < ~, a constant C should 
be determined at other auxiliary conditions In our 



4 



approach it is convenient to redefine the constant C as 
follows C = — A 4 (f — C). Taking into account that p v ~ 
and 4> v oo we transcribe Eq. ( 26 ) into the form 



i(S)'- A '( c - 



A" 



(27) 



Since the second derivative of a chameleon field in the 
spatial region z 2 < ^ is always negative 



n+1 



(28) 



a chameleon field can never reach a minimum between 
two mirrors. This agrees well with the assumption that 
a chameleon field reaches a minimum only at z — > ±oo. 
Then, due to a symmetry of the spatial regions invariant 
under a transformation z < — > — z a chameleon field as a 
scalar field should satisfy the constraint <p(z) — <p(—z). 
Below following [3] we solve a problem of a chameleon 
field between two mirrors assuming that a derivative of 
a chameleon field at z = is continuous. 



A. The chameleon field in the spatial region 
z 2 < with continuous derivative at z = 



Due to a symmetry of a spatial region z 2 < 4r a re- 
quirement of a continuity of a derivative of a chameleon 
field assumes that 



0. 



z=0 



(29) 



we obtain the constant C 
C = 



Setting z = in Eq.(27) and using the constraint Eq.(29) 

A 4+n 



(30) 



we transcribe Eq. ( 32 ) into the form 

rtU) 



n 



z 



dt 



2(n + I)A T J {1 _ t 2 ) ^ 2 



(34) 



where 1/A = A v /n(n + f )(A/0 O )("+ 2 )/ 2 and i(0) 



The integral in the r.h.s. of Eq.(34| can be represented 



in a simpler form. After the change t — sin?/; for < ?/; < 
^ we arrive at the equation 



77, Z i y \rs 2^ 

— — — — =T/ cos* ip dip, (35) 
2(n + 1) A Jo 



where ^((^) 



arcsm 



inv'f - r/^o 1 - 



The integral in the r.h.s. of Eq.(34 1 can be represented 



in terms of the incomplete Beta function [5] . For this aim 
we use a relation 



t(<t>) 



dt 



(1 - t 2 ) 2n 



\B(t 2 (cb) 



I n + 2 
2' 2n 



(36) 



In terms of the incomplete Beta function Eq. ( 34 1 reads 



2(n + I) A 



I n 



2n 



(37) 



From the boundary conditions for a chameleon field at 
z = =p| we obtain a relation between 0^ and (/>o 



2(n + I) A 



SI 



1 n 



2?t 



(38) 



Using the boundary conditions for the derivatives of a 
chameleon field we relate <f>d to 4> m and <pQ as follows 



f 



hn+l 



(39) 



where 0o = <^(0). As a result Eq.(27l takes the form 

A" 



f 



2 (fr=n 



/A r ' 



(31) 



A solution of this equation can be written in the following 
standard form 



± 



1 



n+2 

A~ 



if 2 dip 



(32) 



Suppose that 0o does not depend on /3 in the strong 
coupling limit ft 3> 1. We confirm such a property of <j)Q 
below for the exact solution obtained for n — 2. We also 
show that for the experimental density p m = 2.51 g/cm 3 
of a mirror and d — 25.5 /im [5] the independence of 
of the coupling constant j3 starts for f3 > 10 5 . Since in 
the strong coupling limit (f> m ^> </>^ +1 / <j>Q , <j>d is related 
to d>m as 



(40) 



For the subsequent analysis of the solution of Eq. (32 1 we 



propose to make the following changes of variables. First, 
making a change of variables 



t = Jl- 



(33) 



Substituting Eq. ( 40 1 into Eq. ( 38 ) and taking into account 



that in the strong coupling limit 4>q 3> 4>d we arrive at 
the equation 



2rc 



2(n + l) A 



(41) 



5 



For the derivation of Eq. ( 41 ) we have used a relation |5] 

1\ fn + 2- 



B 1 



1 n + 2 
2' 2n 



2n 



1 



(42) 



where T(a) is Euler's Gamma function and T(l/2) = 
[5]. Solving Eq.(41 ) with respect to 4>o we obtain 4> as a 



function of A, d and n only 



A 



2 

n+2 



V 2n J 



Ad 



(43) 



Thus, we may assert that in the spatial region z 2 < 4- 
in the strong coupling limit /3 > 10 5 , being valid for the 
experimental density of a mirror p rn = 2.51 g/cm 3 [S]and 
d = 25.5 /im , any solution of a chameleon field with a 
continuous derivative at z = should depend on two 
parameters (fid = <j> m and (f> Q , determined by Eq.Q 
for p — p m and Eq. ( 43 ) , respectively. This assertion we 
prove below for the exactly solvable case n = 2. 



B. The chameleon field in the spatial region 
z < ^- with continuous derivative at z = 0. Exact 
solution for n = 2 



32 1 for n = 2 the integral over 



As it is seen from Eq 
if can be calculated in therms of elementary functions. 
Setting n — 2 and integrating over ip we obtain 



V2Az 



This defines 



A 2 

as a function of z 



1 



(*) = 0o(l 



1 rV/2 



3 Ag 



(44) 



(45) 



where 1/A = Av / 6(A/0 O ) 2 . At the next step we have 
to express the parameter </>o in terms of 4> c i = 4>(±d/2). 
Setting z = ±| we obtain 



(46) 



where 1/Ad = Av / 6(A/i/) ( j) 2 . The solution of this alge- 
braical equation with respect to §\ s 





(47) 



Thus, the solution for a chameleon held, defined for n 
2, is 



(z) = <f>, 




(48) 



where A^ is defined by 



Using the relation 




(49) 




one can show that the solution Eq. ( 48 ) satisfies Eq. ( 28 ) 
at n = 2. 

From the boundary conditions for the first derivatives 
of a chameleon field we may define the parameter (fid in 
terms of 6 m 



(51) 



If 0o does not depend on the coupling constant /3 the 
second term in the r.h.s. of Eq.(51 1, taken in the strong 
coupling limit /3 > 10 5 , should be smaller compared with 
the first one. Neglecting the contribution of the second 
term we obtain that in the strong coupling limit (fid can 
be approximated by (fid = \(fi m - This agrees well with 



Eq.(40) at n = 2. Since for n = 2 the problem is exactly 



solvable, we may check such a supposition by using the 
exact relation between (fid and <po, given by Eq.(47), and 



analysing the solution of Eq.(51 numerically. 



For this aim we transcribe Eq.(51 1 into the form 



A3 




(52) 



3A^ 



where we have used Eq.(47|. Denoting X = (fid/4> m we 
rewrite Eq.(52) as follows 



1 



1 + 4/1 



d 2 1 



0, (53) 



where X m is defined by 

A m \<b m ) V 2AI A 3 / 



(54) 



() 



For numerical analysis we use p m = 2.51 g/cm 3 and d = 
25.5 /zm [5]. One can show that the function f(X) has a 
real root, which for f3 > 10 5 is practically equal to X = 
3/2. This confirms our supposition that the parameter 
fa does not depend on the coupling constant /3 in the 
strong coupling limit. 

In the strong coupling limit /3 > 10 5 , where fa — 
(3/2) <f> m , the ratio 4>d/(f>o ^ s given by 



2^3^ 
a 



2V3 1 



d A^6A 2 



Ad A 2 



(55) 



Since (j) 2 , is cancelled, we obtain fa as a function of A and 
d 



5.9 /im, which is a natural scale for neutron quantum 
states in the gravitational field of the Earth [6l [7] . In- 
deed, since a contribution of a chameleon field to the 
gravitational potential of the interaction with ultra-cold 
neutrons to the matrix elements (p\4>\p) for low-lying 
gravitational states \p) — |1) and |3) is localised around 
(| + z) ~ Iq due to the Airy functions [3 [7], the solution 
Eq.(58), obtained in the strong coupling limit f3 > 10 5 , 
should be valid around (| + z) ~ £ . This agrees also 
well with [3J. 



A 

2V4 



'Ad. 



(56) 



This confirms our assertion about the independence of fa 
of the coupling constant /3 in the strong coupling limit 
and, correspondingly, Eq.(43l for n = 2. The scale Ao is 
equal to 



A = 



1 



1 1 



AV& A 2 A-y/6 A 2 ^2 



1 A 2 (Ad) = ' ; 



2-y/3' 



(57) 



Thus, in the strong coupling limit the solution Eq.(48) 
can be transcribed into the form 



A 

2V4 

2 3/4 



Add 

A , ,. 

A" 



d 2 



1/2 



1/2 



(58) 



C. The chameleon field in the spatial region z < ^- 
with continuous derivative at z = 0. Approximate 
solution of non— linear equation for arbitrary n 



Now we may proceed to solving Eq.(|34[) for arbitrary 
n. The only approximation, which may lead to an an- 
alytical solution of the problem, is t <C 1. This allows 
to get an agreement with the exact solution at n = 2. 
Such an approximation gives us the so-cold low-contrast 
solution [2]. Neglecting in Eq.(34| the term t 2 with re- 
spect to a unity and integrating over t we arrive at the 
equation 



y/2Az 



A 



(60) 



In the vicinity of z ~ =f|, where the contribution of the 
mirror, localised at z = ±|, can be neglected and the 
problem under consideration reduces to the problem of a 
chameleon field coupled to ultra-cold neutrons, bouncing 
in the gravitational field of the Earth above a mirror, the 
solution Eq. ( 58 1 takes the form 



0(z) = 2 3 / 4 A(A(; 



±z 



1/2 



(59) 



taken in the 



It agrees well with the solution Eq.(22) 
strong coupling limit. 

As we have shown the exact solution of a chameleon 
field in the spatial region z 2 < ^- , carried out for n = 2, 
confirms fully our suppositions, concerning a dependence 
of the parameters 4>d and 4>$ on the coupling constant /3 
in the strong coupling limit. 

The solution Eq. ( 58 ) can be applied to the analysis of 



the contribution of a chameleon field to the transition 
frequencies of the quantum gravitational states of ultra- 
cold neutrons bouncing in the gravitational field between 
two mirrors. Following [3J, the gravitational potential, 
corrected by the contribution of a chameleon field, takes 



the form Eq.(21 



, ready pointed out that due to the Airy 
wave functions the solution of a chameleon field is lo- 

:| + z) ~ to = (2m 2 gr^ = 



As we have a 



calised in the vicinity of (I 



Solving this equation with respect to </> we obtain a solu- 
tion for a chameleon field 



(z)= 4 i -27^TT)Ay j 



(61) 



where 1/A = Ajn(n + l)(A/(/> )™+ 2 . At n = 2 Eq.(61| 
reduces to Eq.((45l. 



Now we have to replace the parameter (fi by fa. Set- 
ting z = ±| we obtain 



(n + l)A 2 0" 



d 2 ^r 2 

n+2 ' 



(62) 



where 1/Xa = Ayjn(n + 1){A/ 4>d) n+2 ■ For the subse- 
quent calculations it is convenient to rewrite Eq. ( 62 1 as 
follows 



(63) 



1 + 



(n + l)A 2 2 



7 



Using Eq.(63) we transcribe Eq.(61) into the form 



/ n z 2 \ V« 

J) = 0o l 1 - 2 (^TI) AgJ 

n z 2 <$f\h 
2(^+1) 



= 00(^1- 
= <M 1 



l/n 



n d 2 </> 2 



(n + l)A 2 2 



d 2 



x 1 + 



8(n + l)A 2 0g 

n 1 /d 2 
2(n+l)A 2 VT 



M 1 + ^TT)All 




(64) 



where we have denoted = Xd{4>o/ 4>d)- Thus, an ap- 
proximate solution for a chameleon field for an arbitrary 
n is 




For n = 2 the function Eq.(61) reduces to the function 
Eq.gsl 



From Eq. ( 62 ) in the strong coupling limit we define 



by the expression 



A —= Ad 
-2^/2 



2 
n+2 



(66) 



For 1 < n < 10 such an expression reproduces the exact 
result Eq.(43) with an accuracy better than 6.3%. Thus, 
we may assert that the property of <po to be indepen- 
dent of the coupling constant j3 in the strong coupling 
limit is a general property, which does not depend on the 
approximation, but an exact dependence of 4>o on the 
parameters A. d and n depends, of course, on it. 



In order to show that the function Eq.(65) can be 



used as a solution for a chameleon field in the spatial 
region z 2 < ^- and satisfies the equation of motion for a 
chameleon field, we have to derive the equation of motion 



Eq. ( 28 1 at the same approximation, which we have used 



for the derivation of Eq. ( 65 1 . For this aim we have to de- 



rive the equation of motion for t((f>). The first derivative 
of <fr is equal to 



dz 



_ i dt 

dz 



(67) 



and the second one is 
d 2 6 2 



dz 2 



<M(i-0' 



i ft 

n dz 2 
2 , , .I_i/dt\ 2 

-Mi-t 2 ) n 

2 2(n 
n 



ii 



-^o(l-t 2 )« 



2 

\dz/ 



(68) 



The r.h.s. of Eq.(28) can be transcribed into the form 

A n+4 nA «+4 ! 



77 



An+1 



rt+1 



(69) 



Thus, the equation of motion Eq.(28), rewritten for t{<j>), 
takes the form 



t d 2 t / dt\ 2 
dz 2 V dz i 



2(77-1) t 2 



/dt\ 2 



1 - t 2 \dz 



in+l 

b o (1 



(70) 



Neglecting the contributions of the terms of order 0(t 2 ) 
we reduce Eq. ( 70 ) to the form 



d 2 t / dt n 2 
dz 2 V <iz / 



A n+4 
00 



(65) In our approximation the solution of Eq.(34) is 



ra+4 
77 A 2 



(71) 



(72) 



Substituting Eq.(72) into Eq.(72| we satisfy the approx- 



imate equation of motion. This means that we may use 
the solution Eq. ( 65 ) for a chameleon field in the spatial 

region z 2 < 

Since we are interested in the solution for a chameleon 
field in the strong coupling limit, using the rela tion s 4>d — 
4> m and Eq. ( 66 ) we reduce the solution Eq. ( 65 ) to the 
form 



n— 4 



A 2™(™+2) n n+2 (Ad) n (™+2) 

4 



x 



(A-(! 



Z 



(73) 



We would like to remind that the solution Eq.(73) is 
valid for j3 > 10 5 if the mirror density is equal to 
p m = 2.51 g/cm 3 and d = 25.5 g/cm 3 [5]. For n = 2 
it coincides with Eq.( |58[ ). In the vicinity of the mirrors 
z = =F | j the solution Eq. ( 73 ) takes the form 



n— 4 



z) = A2™( n + 2 )77«+2 (Ad)»+2 (k{^- 



±z 



(74) 



One may see that the solution Eq.(74|) reproduces the 
solution for a chameleon field above (1 
Eq.j22l) only for n = 2. 



Delow) a mirror (see 



We would like to note that it is obvious that even if the 
solution for a chameleon field with n = 2 and arbitrary 
coupling constant j3 does not reproduce in the vicinity 
of a mirror the solution for a chameleon field above one 



mirror. This means that fact that the solution Eq.(73) 



does not reproduce in the vicinity of a mirror the solution 
for a chameleon field above a mirror (see Eq.(22| should 



not be evaluated as a strong argument against to use such 
a solution for the problem under consideration. 



That is why in section IV we apply the solution Eq.( 73 ) 



to the estimate of the coupling constant (3 from the con- 
tribution of a chameleon field to the transition frequen- 
cies of quantum gravitational states of ultra-cold neu- 
trons, bouncing in the gravitational field of the Earth 
between two mirrors [5]. 



D. The chameleon held in the spatial region 
z 2 < with continuous derivative at z = 0. Solution 
of the linearised equation of motion for arbitrary n 



In order to complete the analysis of solutions for a 
chameleon field, possessing a continuous derivative at z = 
0, we have to consider the linearised equations of m otion. 
In the spatial region z 2 < ^- the linearised Eq.(28) takes 
the form 



Since we are interesting in the solution for a chameleon 
field in the strong coupling limit, the solution Eq.(79) 



may be valid in the strong coupling limit if it satisfies 
the following constraint. Indeed, in the strong coupling 
limit (pd is proportional to <p m , i.e. <fid — </> m , and 
commen sura ble with zero. This implies that the solu- 



tion Eq.(79| should vanish at z = ±| (see Eq.(58) and 
Eq.(73)). Setting z = ±f we obtain a constraint on (f> 



= A| 



Ad 



2 
n+2 



in 



(80) 



One can show that for 1 < n < 10 the parameter 



defined by Eq.(80), fits the exact expression Eq.(43) with 



an accuracy better than 8.7%. 



The solution Eq.(80l with cj> , defined by Eq.(80), can 



be also applied to the estimate of the coupling constant f3 
from the contribution of a chameleon field to the transi- 
tion frequencies of quantum gravitational states of ultra- 
cold neutrons, bouncing in the gravitational field of the 
Earth between two mirrors [5]. 



d 2 ^ 
dz 2 



A 4+« A 4+n 



(75) 



IV. NUMERICAL SOLUTION OF THE 
PROBLEM 



Denoting c£>q — <f) — <p we define the following linearised 
equation for a chameleon field in the spatial region z 2 < 



d 2 ^ 
dz 2 



\2 
A 



iA 4 +" 

Tn+l ' 



where 1/Ag = n(n + l)A"+ 4 /^ i+2 . As <j>(z) = 
solution of Eq. (l76|) is 



<p(z) 



2 nA 4+n 

J,n+1 
00 



A cosh 



(i). 



(76) 
-z) the 

(77) 



where A is a constant of integration. Since </?(0) = 0, we 
get 



(p(z) 



A 4+n 

\2 

An+l 



1 



x sinh 2 ( - Z 



cosh 
2 



= 9o 



An 
sinh 2 



- A o 



2nA 4+r ' 

rn+l 



Z 



\2X J T "n + l 
This gives a chameleon field <f>(z) equal to 

2 



(«) = 0o(l- 



1 



sinh 



2A 



(78) 



(79) 



This solution agrees well with that obtained in (J (see 
Eq.Q of H at VI = 0). 



I n th is section, we propose a numerical solution of 



Eq. ( 34 ) for a chameleon field in the spatial region z 2 < 



between two mirrors with a density p m = 2.51 g/cm 3 sep- 
arated by the distance d — 25.5 /im [5]. First of all we 
would like note that due to a symmetry of the spatial re- 
gion a chameleon field as a scalar field is symmetric with 
respect to a transformation z < — > —z, i.e. 4>{z) — <fi(—z). 
For a numerical solution we transcribe Eq.( |34[ ) into the 
form 



2(n + l) 



N 

A 



dt 

(1 - t 2 ) 2n 



(81) 



where 1/A = Ay/n(n + l)(A/(/> ) Ii 2 _ , 4>q is given by 
Eq.([43l) and \z\is the absolute value of 2. As a result we 
define a chameleon field as a function of \z\, i.e. cj>(z) = 
f(\z\). Such a non-analytical dependence of a chameleon 
field on a spatial variable does not prevent it to satisfy the 
equation of motion Eq.(28). Indeed, the first derivative 
of a chameleon field wit 



1 respect to z is equal to 



d(f>(z) 
dz 



e(z) 



df(\z\) 
d\z\ : 



(82) 



where e(z) is a sign function, defined by d\z\/dz = e(z) — 
9(z) — 6(—z) and 0(±z) are the Heaviside functions |11) . 
For the second derivative we obtain the following expres- 
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Z[ fim] 



FIG. 1: The profiles of a chameleon field, calculated in the 
strong coupling limit as the solutions of Eq. ( |81[ ) in the spatial 

region z 2 < ^ and n € [1, 10]. 



sion 



dz 2 



2S(z) 



d\z\ 



df(\z\) 



d\z\ 



z=0 



+ s 2 (z) 



d\z\ 2 

dV(M) d 2 f(\z\) 



d\z\ 



d\z\ 



(83) 



where we have used that de(z)/dz = 2S(z) and e 2 (z) — 1 
[llj . Thus, due to a continuous derivative of a chameleon 
field at z = a dependence of a chameleon field on the 
absolute value of z does not introduce an additional terms 



violating the equation of motion Eq.(28) 



For n = 1 the integral Eq.(81) can be calculated ex- 
actly. We obtain 



Aq 



= ty 1 — t 2 + arcsini. 



(84) 



This gives a chameleon field as a function of t(|z|/Ao) 



(*) = 6)(l 



An 



(85) 



In Fig. 1 we show the profiles of a chameleon field, ob- 
tained in the spatial region z 2 < 4- in the strong coupling 
limit by a numerical solution of Eq.(81) for the experi- 
mental mirror density p m = 2.51g/crrr and d — 25.5 /im 
and n € [1, 10]. The lower red line gives the solution for 
n = 1, the green line corresponds to the solution for n = 2 
and so on. In Fig. 2 we give the profiles of a chameleon 
field in the 3D picture. In Fig. 3 wWe compare the nu- 
merical solutions of Eq.(81) with the approximate solu- 



tion Eq. ( 73 1 , which we represent in the more convenient 
form 



1 



2(n + l)\ 2 



d? 
4 



(86) 



where 1/Ao = Ai/n(n + l)(A/</) ) 2i 2 _ but O is given by 
Eq.([66|. The dashed lines describe the approximate 
solutions, whereas the solid line correspond to the exact 
numerical solutions of Eq.(81 1. One can see that for n — 



<f>( z)/0o 




FIG. 2: The profiles of a chameleon field, calculated in the 
strong coupling limit as the solutions of Eq.( 81 1 in the spatial 
region z 2 < ^ and n G [1, 10]. 



2 the dashed and solid lines coincide, since for n = 2 the 
problem admits the analytical solution in the elementary 
functions. Such an agreement confirms the correctness 
all numerical solutions, obtained for n £ [1, 10]. Starting 
with n > 3 the approximate solution is larger compared 
with the exact numerical solutions. 

Following such a behaviour of the approximate solu- 
tions and keeping in mind that the exact solutions of a 
chameleon field taken in the vicinity of a mirror z ~ 



should reproduce the solutions Eq.(22), for the fit of the 



exact numerical solutions of a chameleon field we propose 
the function 



»-*•('-£) 



2 

re +2 



(87) 



V2 d\A 



where <po is equal to 



n + 2 , ,\ n+2 



4V2 



Ad 



(89) 



The results of the fit are shown in Fig. 4. One may see 
that Eq. ( 87 1 fits well the numerical solutions of the non- 



linear equation of motion of a chameleon field and can be 
used for the analytical description of a chameleon field 
in the spatial region z < between two mirrors with 
a density p m = 2.51 g/cm 3 in the strong coupling limit 
f3 > 10 5 and an arbitrary n. For n > 2 the parameter 
(f>o, given by Eq.(89) reproduces the exact 4> 0l defined by 



Eq.(43), with an accuracy better than 3.7%. 
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FIG. 3: The profiles (solid lines) of a chameleon field 
4>(z ) / (f) (0) in the spatial region z 2 < =± as the solutions of 
Eq.(81l, calculated numerically in the strong coupling limit, 



in comparison with the approximate solution Eq.(86 (dashed 
lines) for n 6 [1,10]. 
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FIG. 4: The analytical fit of the numerical solutions of the 
non-linear equations of motion of a chameleon field cj>(z) / </>(0) 
in the spatial region z 2 < ^- , obtained in the strong coupling 
limit for n G [1, 10]. 



V. CONCLUSION 

We have analyzed the solutions for a chameleon field, 
localized between two mirrors in the spatial region z 2 < 
and its influence on the transition frequencies of 
the quantum gravitational states of ultra-cold neutron, 
bouncing between two mirrors in the gravitational field 
of the Earth. 

For the power n = 2 of the potential of a chameleon 
field we have found the exact analytical solutions of non- 



linear equations of motion. We have shown that for the 
experimental density of mirrors p m — 2.51 g/cm 3 and a 
relative distance between mirrors d = 25.5 ^m [5] the 
obtained solutions obey the strong coupling regime with 
j3 > 10 5 . In such a regime a chameleon field does not 
depend on the coupling constant. As a result the contri- 
bution of a chameleon field to the transition frequencies 
of ultra-cold neutrons, bouncing in the gravitational field 
of the Earth between two mirrors, is proportional to the 
coupling constant f3. 
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In the strong coupling limit f3 > 10 5 we have obtained 
the approximate solutions of non-linear and linearized 
equations of motion for a chameleon field and compared 
them with the exact numerical solutions of non-linear 
equations of motion and the analytical fit of the exact 
numerical solutions of non-linear equations of motion. 
The analytical fit of the exact numerical solutions is taken 
in the form reproduci ng i n the vicinity of a mirror at z ~ 
the solutions Eq.(22), coinciding with the solutions 
obtained in [3]. For n > 2 the accuracy of the analytical 
fit of the exact numerical solutions is better than 3.82 %. 

Eq. [87] may be used, for example, to calculate bounds 
on the coupling constant ft by comparing the transition 
frequency uj a b with its theoretical expectation uj a \ co : 



.thco 



^({a\4>(z)\a) 



(b\<Kz)\b)). (90) 



Such resonant transitions between quantum states of a 
neutron in the gravity potential have been measured by 
the qBounce Collaboration [5 . So far a chameleon field 
has been calculated in the infinitely large spatial region 
above a mirror, whereas the ultra-cold neutrons bounce 
in the gravitational field of the Earth between two mirrors 
[5] with a relative distance d. 

Finally we would like to mention that the non-linear 
Eq.(27l with C — admits non-linear solutions with 
a discontinuous derivative of a chameleon field. Such 



solutions are equal to the solutions Eq.(20) with a re- 
placement z — > \z\. The first derivative of these solu- 
tions is proportional to the sign function e(z) [11 . Since 
e 2 (z) = 1 llj, the second derivative satisfies Eq.(28) 
with an additional term — An <Jd S(z), where a Dirac 6- 
function 5(z) appears as a derivative of the sign function 
e'(z) = 2S(z) [TTj . Then, An ad is defined by a jump of 
the first derivative at z = 0, where has a dimension 
of surface density of scalar particle. The factor Air is 
introduced by analogy with electrostatic [12] . An anal- 



ogy between the chameleon field theory in the thin-shell 
regime and the electrostatic has been drawn in [13l [14] . 
A deviation from Eq. ( 28 1 by the term —An adS(z) might 



imply that the Hamilton density of a chameleon field in 
the spatial region 



< ~ should be defined as follows 



H(z) = \ ( d ^) 2 + V{<t>{z)) - 4na5(z)<P(z). (91) 

The parameter a can be unambiguously determined due 
to a self-interaction of a chameleon field. Indeed, the 



Hamilton density Eq.(91) defines the equation of motion 
Eq.(28) with the additional term — AnaS(z). Solving this 



equation for the regions — | < z < and < z < 



d 

2> 

respectively, we arrive at the solution Eq. ( 20 ) with a re- 
placement z — > \z\. This gives a — ad- When the dis- 
tance between two mirrors tends to infinity, the surface 
density id vanishes and we arrive at the solutions for a 
chameleon field above (below) one mirror. 

Of course, we do not stand for the reality of such 
solutions due to a necessity to introduce an additional 
term — Ana 5{z) <fi to the Hamilton density, an influence of 
which on the properties of a chameleon field is not clear. 
The solutions with a discontinuous derivative might be 
accepted as an artifact of the solutions of Eq.(27). 
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